<html>
 <head>
  <meta charset="utf-8"/>
  <meta content="width=device-width, initial-scale=1, maximum-scale=1, user-scalable=no" name="viewport"/>
  <title>
   主题：应如何调用这个函数，得到结果？  | 数螺 | NAUT IDEA
  </title>
  <link href="http://cdn.bootcss.com/bootstrap/3.3.6/css/bootstrap-theme.min.css" rel="stylesheet"/>
  <link href="http://cdn.bootcss.com/bootstrap/3.3.6/css/bootstrap.min.css" rel="stylesheet"/>
  <style type="text/css">
   #xmain img {
                  max-width: 100%;
                  display: block;
                  margin-top: 10px;
                  margin-bottom: 10px;
                }

                #xmain p {
                    line-height:150%;
                    font-size: 16px;
                    margin-top: 20px;
                }

                #xmain h2 {
                    font-size: 24px;
                }

                #xmain h3 {
                    font-size: 20px;
                }

                #xmain h4 {
                    font-size: 18px;
                }


                .header {
	           background-color: #0099ff;
	           color: #ffffff;
	           margin-bottom: 20px;
	        }

	        .header p {
                  margin: 0px;
                  padding: 10px 0;
                  display: inline-block;  
                  vertical-align: middle;
                  font-size: 16px;
               }

               .header a {
                 color: white;
               }

              .header img {
                 height: 25px;
              }
  </style>
  <script src="http://cdn.bootcss.com/jquery/3.0.0/jquery.min.js">
  </script>
  <script src="http://cdn.mathjax.org/mathjax/latest/MathJax.js?config=TeX-AMS-MML_HTMLorMML" type="text/javascript">
   MathJax.Hub.Config({elements: ["bbpress-forums"]});
  </script>
  <script src="http://nautstatic-10007657.file.myqcloud.com/static/css/readability.min.js" type="text/javascript">
  </script>
  <script type="text/javascript">
   $(document).ready(function() {
                 var loc = document.location;
                 var uri = {
                  spec: "http://cos.name/cn/topic/107480/",
                  host: "http://cos.name",
                  prePath: "http://cos.name",
                  scheme: "http",
                  pathBase: "http://cos.name/"
                 };
    
                 var documentClone = document.cloneNode(true);
                 var article = new Readability(uri, documentClone).parse();
     
                 document.getElementById("xmain").innerHTML = article.content;
                });
  </script>
  <!-- 1466447873: Accept with keywords: (title(0.25):主题,论坛,函数,结果, topn(0.133333333333):RSS订阅,软件应用,语言,引号,数据类型,代入,帖子,R语言,作者,会员,用户名,结果,主站,新鲜事,参数,原始数据,函数,老师,环境,问题,密码,总计,字符串,版主,错误,主题,普通,论坛,警告,讨论区).-->
 </head>
 <body class="topic bbpress single single-topic postid-107480 single-author sidebar" onload="">
  <div class="header">
   <div class="container">
    <div class="row">
     <div class="col-xs-6 col-sm-6 text-left">
      <a href="/databee">
       <img src="http://nautidea-10007657.cos.myqcloud.com/logo_white.png"/>
      </a>
      <a href="/databee">
       <p>
        数螺
       </p>
      </a>
     </div>
     <div class="hidden-xs col-sm-6 text-right">
      <p>
       致力于数据科学的推广和知识传播
      </p>
     </div>
    </div>
   </div>
  </div>
  <div class="container text-center">
   <h1>
    主题：应如何调用这个函数，得到结果？
   </h1>
  </div>
  <div class="container" id="xmain">
   <div class="hfeed site" id="page">
    <header class="site-header" id="masthead" role="banner">
     <div id="cos-logo">
      <a href="http://cos.name/cn">
       <img src="http://cos.name/cn/wp-content/themes/COS-forest/images/headers/cos-logo.png"/>
      </a>
     </div>
     <div class="navbar" id="navbar">
      <nav class="navigation main-navigation" id="site-navigation" role="navigation">
       <h3 class="menu-toggle">
        菜单
       </h3>
       <div class="menu-%e8%8f%9c%e5%8d%951-container">
        <ul class="nav-menu" id="menu-%e8%8f%9c%e5%8d%951">
         <li class="menu-item menu-item-type-custom menu-item-object-custom menu-item-home menu-item-407772" id="menu-item-407772">
          <a href="http://cos.name/cn/">
           论坛首页
          </a>
         </li>
         <li class="menu-item menu-item-type-custom menu-item-object-custom menu-item-407773" id="menu-item-407773">
          <a href="http://cos.name/cn/forums/">
           讨论区
          </a>
         </li>
         <li class="menu-item menu-item-type-custom menu-item-object-custom menu-item-407774" id="menu-item-407774">
          <a href="http://cos.name/cn/wp-login.php?action=register">
           注册
          </a>
         </li>
         <li class="menu-item menu-item-type-custom menu-item-object-custom menu-item-407819" id="menu-item-407819">
          <a href="http://cos.name/">
           主站
          </a>
         </li>
        </ul>
       </div>
      </nav>
      <!-- #site-navigation -->
     </div>
     <!-- #navbar -->
    </header>
    <!-- #masthead -->
    <div class="site-main" id="main">
     <div class="content-area" id="primary">
      <div class="site-content" id="content" role="main">
       <article class="post-107480 topic type-topic status-publish hentry" id="post-107480">
        <header class="entry-header">
         <h1 class="entry-title">
          应如何调用这个函数，得到结果？
         </h1>
        </header>
        <!-- .entry-header -->
        <div class="entry-content">
         <div id="bbpress-forums">
          <div class="bbp-breadcrumb">
           <p>
            <a class="bbp-breadcrumb-home" href="http://cos.name/cn/">
             COS论坛 | 统计之都
            </a>
            <span class="bbp-breadcrumb-sep">
             ›
            </span>
            <a class="bbp-breadcrumb-root" href="http://cos.name/cn/forums/">
             讨论区
            </a>
            <span class="bbp-breadcrumb-sep">
             ›
            </span>
            <a class="bbp-breadcrumb-forum" href="http://cos.name/cn/forum/software/">
             软件应用
            </a>
            <span class="bbp-breadcrumb-sep">
             ›
            </span>
            <a class="bbp-breadcrumb-forum" href="http://cos.name/cn/forum/software/r-language/">
             S-Plus &amp; R语言
            </a>
            <span class="bbp-breadcrumb-sep">
             ›
            </span>
            <span class="bbp-breadcrumb-current">
             应如何调用这个函数，得到结果？
            </span>
           </p>
          </div>
          <div class="bbp-template-notice info">
           <p class="bbp-topic-description">
            该主题包含 9 条回复，2个帖子，最后由
            <a class="bbp-author-avatar" href="http://cos.name/cn/profile/109572/" rel="nofollow" title="查看ywh19860616的档案">
             <img src="http://sdn.geekzu.org/avatar/8599115e910c111374c57034b9f33bda?s=14&amp;d=monsterid&amp;r=g"/>
            </a>
            <a class="bbp-author-name" href="http://cos.name/cn/profile/109572/" rel="nofollow" title="查看ywh19860616的档案">
             ywh19860616
            </a>
            在
            <a href="http://cos.name/cn/topic/107480/#post-334618" title="回复：应如何调用这个函数，得到结果？">
             3 年, 11 月 之前
            </a>
            更新。
           </p>
          </div>
          <div class="bbp-pagination">
           <div class="bbp-pagination-count">
            查看 10 个帖子 - 1 到 10（总计 10 个）
           </div>
           <div class="bbp-pagination-links">
           </div>
          </div>
          <ul class="forums bbp-replies" id="topic-107480-replies">
           <li class="bbp-header">
            <div class="bbp-reply-author">
             作者
            </div>
            <!-- .bbp-reply-author -->
            <div class="bbp-reply-content">
             帖子
            </div>
            <!-- .bbp-reply-content -->
           </li>
           <!-- .bbp-header -->
           <li class="bbp-body">
            <div class="bbp-reply-header" id="post-107480">
             <div class="bbp-meta">
              <span class="bbp-reply-post-date">
               2012年7月13日 上午11:48
              </span>
              <a class="bbp-reply-permalink" href="http://cos.name/cn/topic/107480/#post-107480">
               1 楼
              </a>
              <span class="bbp-admin-links">
              </span>
             </div>
             <!-- .bbp-meta -->
            </div>
            <!-- #post-107480 -->
            <div class="odd bbp-parent-forum-999 bbp-parent-topic-107480 bbp-reply-position-1 user-id-109572 topic-author post-107480 topic type-topic status-publish hentry">
             <div class="bbp-reply-author">
              <a class="bbp-author-avatar" href="http://cos.name/cn/profile/109572/" rel="nofollow" title="查看ywh19860616的档案">
               <img src="http://sdn.geekzu.org/avatar/8599115e910c111374c57034b9f33bda?s=80&amp;d=monsterid&amp;r=g"/>
              </a>
              <br/>
              <a class="bbp-author-name" href="http://cos.name/cn/profile/109572/" rel="nofollow" title="查看ywh19860616的档案">
               ywh19860616
              </a>
              <br/>
              <div class="bbp-author-role">
               普通会员
              </div>
             </div>
             <!-- .bbp-reply-author -->
             <div class="bbp-reply-content">
              <p>
               各位前辈好，请教一个问题：
               <br/>
               以下我有一个函数fit701
              </p>
              <p>
               fit701=function(xv,yv,etx,dtx,wa){
               <br/>
               # Model M1
               <br/>
               # Lee-Carter model
               <br/>
               # log m(t,x) = beta1(x) + beta2(x).kappa2(t) + Poisson error
              </p>
              <p>
               # Inputs:
               <br/>
               #   xv = vector of ages, length n
               <br/>
               #   yv = vector of years, length m
               <br/>
               #   etx = m x n matrix of exposures
               <br/>
               #   dtx = m x n matrix of deaths
               <br/>
               #   wa = m x n matrix of weights (0 or 1)
              </p>
              <p>
               mtx=dtx/etx	  # matrix of death rates
               <br/>
               qtx=1-exp(-mtx) # matrix of mortality rates
              </p>
              <p>
               if(max(xv) &gt; 89)
               <br/>
               {
               <br/>
               cat(“Upper age too high – suggest abort programme\n”)
               <br/>
               }
              </p>
              <p>
               n=length(xv)	# number of ages
               <br/>
               m=length(yv)	# number of years
              </p>
              <p>
               cy=(yv[1]-xv[n]):(yv[m]-xv[1])  # cohort approximate years of birth
              </p>
              <p>
               # initialise parameter vectors
               <br/>
               beta1v=(1:n)*0
               <br/>
               beta2v=(1:n)*0
               <br/>
               beta3v=(1:n)*0		# dummy vector, this will stay at 0
               <br/>
               kappa2v=(1:m)*0
               <br/>
               gamma3v=(1:(n+m-1))*0	# dummy vector, this will stay at 0
               <br/>
               ia=array((1:m),c(m,n))	# matrix of year indexes, i, for the data
               <br/>
               ja=t(array((1:n),c(n,m)))	# matrix of age indexes, j, for the data
               <br/>
               ya=ia-ja		 	# matrix of year of birth indexes for the data
               <br/>
               imj=(1-n):(m-1)		# the range of values taken by i-j
               <br/>
               lg=n+m-1		 	# number of different values taken by i-j
               <br/>
               ca=ya+yv[1]-xv[1]		# matrix of years of birth
              </p>
              <p>
               # Now set weights to zero for cohorts with fewer than 5 observations
               <br/>
               for(k in 1:lg)
               <br/>
               {
               <br/>
               nk=sum((ca == cy[k])*wa)
               <br/>
               if(nk &lt; 5)
               <br/>
               {
               <br/>
               wa=wa*(1- (ca == cy[k]))
               <br/>
               }
               <br/>
               }
              </p>
              <p>
               ww=cy*0+1	 # this is a vector of 1’s and 0’s with
               <br/>
               # a 0 if the cohort is completely excluded
               <br/>
               for(k in 1:lg)
               <br/>
               {
               <br/>
               ww[k]=ww[k]*(sum((ca == cy[k])*wa) &gt; 0)
               <br/>
               }
              </p>
              <p>
               # Stage 0
               <br/>
               # Gives initial estimates for beta1(x), beta2(x) and kappa2(t)
               <br/>
               mx=mean(xv)
               <br/>
               for(j in 1:n)
               <br/>
               {
               <br/>
               beta1v[j]=sum(log(mtx[,j])*wa[,j])/sum(wa[,j])
               <br/>
               beta2v[j]=1/n
               <br/>
               }
               <br/>
               kappa2v=(m:1)-(m+1)/2
              </p>
              <p>
               # Stage 1: iterate
               <br/>
               l0=-1000000
               <br/>
               l1=-999999
               <br/>
               iteration=0
               <br/>
               # l1 is the latest estimate of the log-likelihood
               <br/>
               # l0 is the previous estimate
               <br/>
               # we continue to iterate if the improvement in log-likelihood
               <br/>
               # exceeds 0.0001
               <br/>
               while(abs(l1-l0) &gt; 0.0001)
               <br/>
               {
               <br/>
               iteration=iteration+1
              </p>
              <p>
               l0=l1
               <br/>
               # Stage 1B optimise over the beta2(x)
               <br/>
               for(j in 1:n)
               <br/>
               {
               <br/>
               # cycle through the range of years
               <br/>
               dv=dtx[,j]	# actual deaths
               <br/>
               ev=etx[,j]	# exposure
               <br/>
               beta2v[j]=llmaxM2B(beta1v[j],beta2v[j],beta3v[j],
               <br/>
               kappa2v,gamma3v[(n+1-j):(n+m-j)],dv,ev,wv=wa[,j])
               <br/>
               }
              </p>
              <p>
               mhat=mtx*0
               <br/>
               for(i in 1:m)
               <br/>
               {
               <br/>
               mhat
               <em class="d4pbbc-italic">
               </em>
               =exp(beta1v+beta2v*kappa2v
               <em class="d4pbbc-italic">
               </em>
               +beta3v*gamma3v[(n+i-1):i])
               <br/>
               }
               <br/>
               epsilon=(dtx-etx*mhat)/sqrt(etx*mhat)
               <br/>
               l1=sum((dtx*log(etx*mhat)-etx*mhat-lgamma(dtx+1))*wa)
               <br/>
               cat(l1,”-&gt; “)
              </p>
              <p>
               # Stage 1D optimise over the kappa2(t)
               <br/>
               for(i in 1:m)
               <br/>
               {
               <br/>
               # cycle through the range of years
               <br/>
               dv=dtx
               <em class="d4pbbc-italic">
               </em>
               # actual deaths
               <br/>
               ev=etx
               <em class="d4pbbc-italic">
               </em>
               # exposure
               <br/>
               kappa2v
               <em class="d4pbbc-italic">
               </em>
               =llmaxM2D(beta1v,beta2v,beta3v,
               <br/>
               kappa2v
               <em class="d4pbbc-italic">
               </em>
               ,gamma3v[(n+i-1):i],dv,ev,wv=wa
               <em class="d4pbbc-italic">
               </em>
               )
               <br/>
               }
              </p>
              <p>
               # Now apply the constraints
               <br/>
               fac21=mean(kappa2v)
               <br/>
               fac22=sum(beta2v)
               <br/>
               kappa2v=fac22*(kappa2v-fac21)    # ensures that the kappas sum to 0
               <br/>
               beta2v=beta2v/fac22		     # ensures that the beta2’s sum to 1
               <br/>
               beta1v=beta1v+beta2v*fac22*fac21 # =&gt; beta1 needs to be adjusted to compensate
              </p>
              <p>
               mhat=mtx*0
               <br/>
               for(i in 1:m)
               <br/>
               {
               <br/>
               mhat
               <em class="d4pbbc-italic">
               </em>
               =exp(beta1v+beta2v*kappa2v
               <em class="d4pbbc-italic">
               </em>
               +beta3v*gamma3v[(n+i-1):i])
               <br/>
               }
               <br/>
               epsilon=(dtx-etx*mhat)/sqrt(etx*mhat)
               <br/>
               l1=sum((dtx*log(etx*mhat)-etx*mhat-lgamma(dtx+1))*wa)
               <br/>
               cat(l1,” -&gt;”)
              </p>
              <p>
               # Stage 1A optimise over the beta1(x)
               <br/>
               for(j in 1:n)
               <br/>
               {
               <br/>
               # cycle through the range of years
               <br/>
               wv=1	    # can be set to a vector of weights
               <br/>
               # to e.g. exclude duff years
               <br/>
               wv=wa[,j]
               <br/>
               s1=sum(wv*dtx[,j])
               <br/>
               s2=sum(wv*etx[,j]*exp(beta2v[j]*kappa2v+beta3v[j]*gamma3v[(n+1-j):(n+m-j)]))
               <br/>
               beta1v[j]=log(s1)-log(s2)
               <br/>
               }
              </p>
              <p>
               mhat=mtx*0
               <br/>
               for(i in 1:m)
               <br/>
               {
               <br/>
               mhat
               <em class="d4pbbc-italic">
               </em>
               =exp(beta1v+beta2v*kappa2v
               <em class="d4pbbc-italic">
               </em>
               +beta3v*gamma3v[(n+i-1):i])
               <br/>
               }
               <br/>
               epsilon=(dtx-etx*mhat)/sqrt(etx*mhat)
               <br/>
               l1=sum((dtx*log(etx*mhat)-etx*mhat-lgamma(dtx+1))*wa)
               <br/>
               cat(l1,”\n”)
              </p>
              <p>
               }		 # end while loop
              </p>
              <p>
               # calculate number of parameters and deduct 4 for the number of constraints
               <br/>
               npar=length(beta1v)+length(beta2v)+length(kappa2v)-2
              </p>
              <p>
               # Calculate the BIC
               <br/>
               BIC=l1-0.5*log(sum(wa))*npar
              </p>
              <p>
               list(beta1=beta1v,beta2=beta2v,beta3=beta3v,
               <br/>
               kappa2=kappa2v,gamma3=gamma3v,x=xv,y=yv,cy=cy,
               <br/>
               wa=wa,epsilon=epsilon,mhat=mhat,ll=l1,BIC=BIC,npar=npar,mtxLastYear = mtx[m,])
              </p>
              <p>
               }
               <br/>
               plot(beta2v)
              </p>
              <p>
               我的调用方式如下：
               <br/>
               x = as.integer(c(60, 61, 62, 63, 64, 65, 66, 67,
               <br/>
               68, 69, 70, 71, 72, 73, 74, 75, 76, 77, 78, 79, 80, 81, 82, 83,
               <br/>
               84, 85, 86, 87, 88, 89))
               <br/>
               y = as.integer(c(1968, 1969, 1970,
               <br/>
               1971, 1972, 1973, 1974, 1975, 1976, 1977, 1978, 1979, 1980, 1981,
               <br/>
               1982, 1983, 1984, 1985, 1986, 1987, 1988, 1989, 1990, 1991, 1992,
               <br/>
               1993, 1994, 1995, 1996, 1997, 1998, 1999, 2000, 2001, 2002, 2003
               <br/>
               ))
               <br/>
               etx = structure(c(856167, 866287, 903314, 875491, 880448,
               <br/>
               892288, 897373, 906129, 891513, 886044, 905800, 940496, 1021278,
               <br/>
               1062928, 1075375, 1037117, 1056018, 1052845, 1020308, 1024275,
               <br/>
               1035615, 1012013, 1011827, 994493, 949037, 939460, 933858, 995398,
               <br/>
               978947, 984941, 1037744, 1054359, 1119331, 1138132, 1215153,
               <br/>
               1367422, 813957, 832384, 841000, 886880, 865567, 876908, 896311,
               <br/>
               909091, 925398, 916932, 918861, 947455, 969148, 994121, 1035118,
               <br/>
               1047506, 1009979, 1028611, 1025401, 993776, 997940, 1009834,
               <br/>
               984232, 993307, 976781, 932267, 922196, 917577, 977659, 962298,
               <br/>
               968349, 1020164, 1066360, 1105452, 1123898, 1199917, 769980,
               <br/>
               793648, 809970, 825184, 871486, 850010, 861881, 882401, 895572,
               <br/>
               912479, 904679, 907038, 927024, 920673, 947893, 991340, 1007339,
               <br/>
               975370, 997039, 998294, 971480, 980136, 999158, 971062, 980626,
               <br/>
               964923, 920409, 911400, 906473, 966663, 952002, 958203, 1041131,
               <br/>
               1051307, 1090074, 1108149, 730887, 754131, 761967, 797520, 812424,
               <br/>
               859139, 837780, 850030, 870859, 884269, 901900, 893952, 900452,
               <br/>
               928482, 922297, 949057, 992613, 1008397, 976113, 997755, 998713,
               <br/>
               971933, 980666, 1003543, 975933, 986025, 969719, 925486, 915683,
               <br/>
               911581, 972711, 957676, 978987, 1025291, 1035410, 1073948, 694120,
               <br/>
               716735, 735620, 748868, 783125, 796414, 842754, 820393, 831159,
               <br/>
               851027, 863458, 880029, 879452, 900385, 928190, 922262, 948390,
               <br/>
               991946, 1007581, 975680, 997198, 998365, 971835, 981866, 1006136,
               <br/>
               979040, 989062, 973764, 928910, 919807, 916406, 977380, 957302,
               <br/>
               962735, 1008859, 1018693, 660800, 678637, 706601, 720541, 732391,
               <br/>
               765768, 778637, 824789, 801253, 811707, 831060, 843141, 862694,
               <br/>
               864906, 883432, 907341, 898438, 920579, 959249, 971271, 937508,
               <br/>
               955284, 959793, 939209, 949592, 974070, 947836, 958501, 943473,
               <br/>
               900478, 892028, 888601, 953199, 939883, 945434, 991349, 633685,
               <br/>
               643455, 666122, 688795, 701209, 711768, 744987, 757324, 803146,
               <br/>
               779215, 788633, 807738, 816708, 822067, 826637, 846701, 870571,
               <br/>
               864505, 886836, 926358, 939387, 908921, 928135, 932209, 912993,
               <br/>
               923548, 947407, 922890, 933141, 919059, 877065, 868563, 872243,
               <br/>
               934290, 921171, 927025, 610486, 614039, 629797, 645396, 666993,
               <br/>
               678635, 689230, 722362, 734266, 780124, 755703, 765708, 787301,
               <br/>
               796780, 803904, 809868, 831355, 856847, 852011, 875706, 916458,
               <br/>
               931259, 897058, 909180, 913663, 894941, 905034, 929745, 905727,
               <br/>
               916601, 902826, 861400, 868256, 853229, 914644, 901758, 583782,
               <br/>
               586915, 570692, 605096, 619979, 641386, 653640, 664590, 698043,
               <br/>
               710755, 757283, 733811, 743285, 752597, 759744, 763839, 766951,
               <br/>
               786160, 807024, 800260, 819592, 855612, 872167, 847508, 859081,
               <br/>
               863602, 845901, 856738, 880054, 858040, 868975, 855817, 856238,
               <br/>
               847663, 832930, 893913, 553129, 557079, 565585, 544330, 577926,
               <br/>
               592750, 614323, 627086, 638650, 672465, 685353, 731561, 705800,
               <br/>
               719835, 731507, 740038, 745870, 752289, 773836, 796830, 792965,
               <br/>
               814900, 855420, 862899, 838560, 850050, 854363, 837706, 848156,
               <br/>
               872633, 851424, 862324, 844657, 833686, 825860, 811629, 521770,
               <br/>
               524742, 546846, 531810, 505837, 533771, 543806, 559832, 567126,
               <br/>
               573232, 599359, 606474, 657413, 670764, 684652, 695130, 702755,
               <br/>
               708533, 714563, 736428, 757770, 754756, 777624, 829298, 836574,
               <br/>
               813103, 823991, 829034, 812563, 823632, 848179, 827664, 847845,
               <br/>
               819982, 810102, 802793, 493069, 493613, 489228, 521292, 509924,
               <br/>
               487698, 521330, 536585, 558895, 571906, 582857, 617382, 615889,
               <br/>
               629391, 643271, 656880, 667056, 674891, 680522, 687546, 709690,
               <br/>
               732227, 727327, 751182, 801706, 809020, 785910, 797618, 802053,
               <br/>
               786636, 798045, 821777, 798523, 821496, 794293, 785218, 464965,
               <br/>
               462899, 449521, 460001, 489586, 476858, 454388, 486414, 500451,
               <br/>
               520873, 531577, 541572, 579817, 590133, 601112, 611638, 621567,
               <br/>
               628331, 632384, 635025, 638853, 658325, 688165, 711376, 734762,
               <br/>
               784702, 791570, 770146, 781352, 786237, 771921, 783023, 786504,
               <br/>
               770509, 794034, 767594, 435941, 432637, 428927, 421251, 432833,
               <br/>
               464833, 454759, 433625, 468772, 485424, 508535, 522305, 524287,
               <br/>
               536570, 549664, 562947, 575846, 588445, 597690, 604940, 610374,
               <br/>
               618475, 629835, 638876, 661118, 682842, 728947, 736425, 716127,
               <br/>
               727165, 732611, 719175, 756705, 756305, 741418, 765096, 408692,
               <br/>
               402218, 407575, 398877, 390529, 402041, 434255, 424518, 404333,
               <br/>
               439415, 455364, 478352, 494175, 506181, 518832, 531733, 544936,
               <br/>
               558247, 570767, 580312, 587293, 594013, 600377, 608298, 617027,
               <br/>
               639073, 659894, 705499, 712462, 693485, 705067, 710352, 715681,
               <br/>
               725129, 725091, 711218, 380807, 375963, 380169, 378506, 369129,
               <br/>
               360766, 372938, 405468, 395893, 376042, 410407, 426790, 450828,
               <br/>
               463599, 476265, 488639, 501616, 515586, 529234, 542437, 551994,
               <br/>
               560445, 570224, 577377, 585864, 594420, 615387, 636575, 680315,
               <br/>
               687559, 670185, 681390, 695365, 682826, 692940, 693424, 348190,
               <br/>
               348880, 345547, 353216, 352107, 343718, 337013, 350924, 386492,
               <br/>
               378492, 356556, 393999, 407457, 419043, 431096, 442350, 452940,
               <br/>
               465081, 477747, 490662, 502117, 511231, 521203, 540528, 547729,
               <br/>
               556190, 564112, 585037, 604955, 647256, 655104, 638556, 651132,
               <br/>
               660748, 649399, 659990, 312353, 317343, 320291, 319364, 326669,
               <br/>
               325651, 318197, 312252, 326384, 362430, 354026, 333083, 368217,
               <br/>
               378617, 390303, 401778, 412498, 422828, 434639, 447633, 460160,
               <br/>
               472377, 481592, 491871, 510544, 517581, 525378, 534002, 553631,
               <br/>
               573151, 614325, 621876, 606071, 615676, 625752, 615365, 278317,
               <br/>
               282976, 281045, 297607, 297834, 307266, 308946, 303543, 298591,
               <br/>
               315332, 356824, 350445, 333174, 340792, 350879, 361371, 371677,
               <br/>
               381554, 390726, 402038, 413924, 426669, 441344, 462018, 472525,
               <br/>
               490868, 497465, 506050, 514196, 533545, 553352, 593222, 573575,
               <br/>
               570641, 580495, 590742, 246589, 250690, 245642, 257870, 272812,
               <br/>
               271793, 280268, 281212, 274785, 269800, 285253, 321181, 302394,
               <br/>
               308751, 315784, 324947, 334281, 343887, 353005, 361697, 371991,
               <br/>
               384214, 397314, 413871, 433410, 443319, 460297, 467521, 475505,
               <br/>
               483671, 503005, 521959, 528796, 537217, 535475, 545338, 218385,
               <br/>
               222159, 235718, 221534, 229560, 240281, 236160, 241397, 238536,
               <br/>
               229347, 221123, 232063, 262589, 268088, 274795, 280870, 289391,
               <br/>
               298159, 307301, 316436, 324212, 335173, 344376, 352360, 367467,
               <br/>
               384802, 393710, 409973, 416410, 423822, 431991, 449167, 476256,
               <br/>
               491351, 500206, 499428, 191831, 194665, 208358, 215780, 203574,
               <br/>
               214408, 229348, 228547, 237745, 238281, 230671, 224253, 225823,
               <br/>
               231064, 237512, 244519, 250579, 259603, 268728, 278865, 288280,
               <br/>
               297598, 307532, 315275, 322913, 336912, 352664, 361854, 376721,
               <br/>
               383001, 390718, 398318, 404709, 440534, 454061, 463344, 167849,
               <br/>
               170176, 170352, 188680, 194666, 180977, 191234, 205654, 203857,
               <br/>
               211889, 211198, 203570, 198428, 203255, 208777, 214984, 221582,
               <br/>
               227259, 236022, 245372, 254897, 265010, 273384, 279758, 286839,
               <br/>
               294185, 306825, 322103, 330490, 344466, 351051, 358174, 373147,
               <br/>
               369929, 405690, 417949, 145946, 148393, 143516, 151984, 169036,
               <br/>
               173413, 158822, 168511, 181773, 179350, 185783, 184670, 180549,
               <br/>
               183544, 188156, 192055, 196646, 201516, 205290, 212744, 219811,
               <br/>
               228050, 237226, 249171, 255501, 262325, 269022, 281504, 295494,
               <br/>
               303525, 317298, 323430, 319547, 340044, 336089, 371368, 125564,
               <br/>
               128851, 125396, 126958, 134037, 149410, 152316, 137075, 145796,
               <br/>
               157832, 154337, 159741, 158816, 161550, 165031, 169425, 172718,
               <br/>
               176646, 180713, 184009, 190366, 197346, 203495, 211185, 221884,
               <br/>
               227864, 233974, 240702, 251805, 264712, 272624, 285042, 280589,
               <br/>
               288454, 308106, 303998, 118294.174369962, 117593.628507295, 119786.337057443,
               <br/>
               121733.854555658, 118728.512804816, 129075.575196409, 131373.365625957,
               <br/>
               124444.967044179, 123933.568564432, 127128.057698194, 132760.446434037,
               <br/>
               129103.597030915, 137321, 134699, 138443, 142392, 146985, 150549,
               <br/>
               154402, 158689, 161412, 168523, 173121, 171462, 177988, 187182,
               <br/>
               192148, 197978, 203740, 213411, 225037, 231900, 247227, 250442,
               <br/>
               258129, 276966, 96295.436345537, 97930.8887291229, 99599.0501603805,
               <br/>
               103013.874737308, 104289.527596505, 102706.409689194, 107279.134553700,
               <br/>
               106880.084172105, 109137.008461453, 105218.464550382, 111001.424178741,
               <br/>
               112990.134277182, 118289, 117771, 115720, 118882, 122214, 126414,
               <br/>
               129315, 132401, 135164, 137533, 144002, 151735, 150320, 156241,
               <br/>
               164259, 169283, 174439, 179759, 188933, 199339, 209266, 218635,
               <br/>
               221409, 228887, 75689.0029792667, 76349.2036237612, 77135.4425731136,
               <br/>
               82459.0604973551, 83557.3942968323, 85712.0491275005, 81768.8507326564,
               <br/>
               84997.8320666383, 90651.550313127, 88742.970268134, 88845.0012768286,
               <br/>
               91053.6725238645, 98712, 100443, 100221, 98062, 100740, 103982,
               <br/>
               108035, 110707, 112906, 115577, 118031, 124972, 131717, 130542,
               <br/>
               135678, 143252, 147619, 152359, 157609, 165712, 174659, 183063,
               <br/>
               191560, 194338, 61829.1830637757, 61437.1759305036, 61920.0832685924,
               <br/>
               65397.0161028321, 67061.6261035384, 70049.9703368882, 70788.534501024,
               <br/>
               66232.1617345858, 72549.7259693481, 74617.7056289286, 75856.2209195565,
               <br/>
               72987.1832050286, 80441, 82144, 83849, 82960, 80274, 82298, 85108,
               <br/>
               88823, 90412, 92164, 93964, 101263, 107237, 113121, 112076, 117038,
               <br/>
               123569, 127522, 132180, 136825, 142566, 150887, 158556, 166688,
               <br/>
               47640.0101531432, 50176.0216600389, 50012.7591166765, 53027.6740841019,
               <br/>
               54877.9829088755, 57648.0040612573, 57745.9615872747, 56352.7878839157,
               <br/>
               58186.7704543532, 59416.6816143498, 64940.397664777, 62562.206616465,
               <br/>
               64320, 66906, 69284, 71351, 70635, 68463, 71073, 74925, 78939,
               <br/>
               81234, 82499, 79650, 85828, 90988, 95962, 95543, 99802, 105534,
               <br/>
               109371, 113474, 114786, 121844, 128943, 136238), .Dim = as.integer(c(36,
               <br/>
               30)))
              </p>
              <p>
               dtx = structure(c(21137, 20362, 21053, 20598, 20766, 21026,
               <br/>
               20378, 19689, 19332, 18437, 18790, 18119, 19212, 19062, 18826,
               <br/>
               18446, 18507, 18680, 17681, 17509, 17440, 16210, 15930, 15402,
               <br/>
               14878, 14538, 14379, 14708, 14141, 13954, 13893, 14016, 14152,
               <br/>
               14542, 15373, 16717, 21093, 21576, 20997, 21426, 22068, 21349,
               <br/>
               21138, 20328, 20488, 19787, 19608, 19382, 19207, 19856, 20094,
               <br/>
               19723, 19317, 19526, 19084, 18710, 18359, 17818, 16884, 16950,
               <br/>
               15977, 15912, 15208, 15199, 15354, 14830, 14780, 14881, 14657,
               <br/>
               14909, 15326, 16369, 22667, 22628, 22831, 22348, 23448, 23410,
               <br/>
               22314, 22219, 22310, 21906, 21298, 20757, 21006, 20407, 21555,
               <br/>
               21573, 21504, 21027, 21023, 20583, 19930, 19423, 19275, 18131,
               <br/>
               17868, 17670, 17082, 16338, 16482, 16170, 15974, 15506, 15832,
               <br/>
               15801, 16305, 16600, 23364, 22818, 23381, 23675, 23894, 24137,
               <br/>
               23362, 22804, 23007, 22916, 22873, 21796, 21637, 21943, 20956,
               <br/>
               22485, 22742, 22476, 21750, 21785, 21673, 20919, 20442, 19971,
               <br/>
               19180, 19154, 18385, 17748, 17441, 16617, 17078, 16626, 16511,
               <br/>
               16672, 16861, 16901, 23909, 23759, 23940, 23855, 25032, 24501,
               <br/>
               24242, 23920, 23576, 23459, 23996, 23420, 23125, 22804, 22615,
               <br/>
               22451, 23349, 23779, 23595, 22722, 22772, 22258, 21588, 21640,
               <br/>
               20894, 20503, 20443, 19449, 18778, 18121, 17514, 18254, 17744,
               <br/>
               17233, 17668, 17686, 25155, 24952, 25534, 25315, 25402, 26151,
               <br/>
               25495, 25265, 25308, 24760, 24643, 24737, 24891, 24157, 23878,
               <br/>
               24163, 23624, 24991, 25225, 24798, 24225, 23691, 23626, 22657,
               <br/>
               22407, 22503, 21216, 20935, 20201, 19320, 18908, 18538, 18614,
               <br/>
               18366, 18106, 18427, 23831, 24553, 24964, 25200, 25584, 25850,
               <br/>
               25891, 25057, 25424, 25273, 25329, 24966, 25689, 25347, 24881,
               <br/>
               24465, 24992, 24617, 25119, 25628, 25696, 24386, 24367, 24049,
               <br/>
               23170, 23352, 23163, 22152, 21759, 21043, 20462, 19743, 18960,
               <br/>
               19340, 18993, 18805, 27003, 24676, 25690, 25921, 26864, 26654,
               <br/>
               26479, 26549, 26216, 26176, 26621, 25567, 26743, 26416, 26337,
               <br/>
               26011, 25668, 26082, 25560, 26453, 26724, 26240, 25331, 25501,
               <br/>
               25109, 24941, 24594, 24254, 23058, 23098, 22022, 21115, 20658,
               <br/>
               19994, 20463, 20079, 27021, 26871, 25138, 25785, 26792, 26957,
               <br/>
               26819, 26882, 26967, 26236, 27039, 26615, 26988, 26943, 27258,
               <br/>
               27080, 26638, 26587, 26986, 26482, 27667, 27512, 26971, 25991,
               <br/>
               26531, 26854, 25485, 25253, 25091, 24331, 23777, 23073, 22082,
               <br/>
               21299, 20769, 21144, 25507, 27135, 27237, 25408, 26978, 27459,
               <br/>
               27157, 26608, 27436, 27428, 27163, 27181, 28148, 27404, 27783,
               <br/>
               28498, 28289, 27988, 27723, 27761, 27619, 28018, 28061, 28183,
               <br/>
               27011, 27795, 27550, 26821, 26551, 25998, 25024, 25060, 23621,
               <br/>
               23111, 22259, 21860, 27167, 26198, 27741, 27687, 26400, 27716,
               <br/>
               27258, 26867, 27577, 27489, 28053, 27411, 28736, 28704, 28172,
               <br/>
               29089, 28964, 29432, 28934, 28155, 28648, 27956, 28751, 29132,
               <br/>
               28810, 28816, 28703, 28552, 28017, 27833, 26962, 25849, 25690,
               <br/>
               24481, 24091, 23025, 27269, 26523, 25626, 27679, 28776, 26401,
               <br/>
               26874, 27039, 26966, 27466, 27921, 27959, 28297, 28977, 29096,
               <br/>
               29153, 29592, 30090, 30102, 29486, 29148, 29060, 28329, 29878,
               <br/>
               30110, 30663, 29723, 29800, 29578, 28528, 28667, 28204, 27007,
               <br/>
               26412, 25839, 24264, 29853, 27868, 27050, 26600, 29170, 29366,
               <br/>
               27147, 27534, 27931, 27909, 28493, 28336, 29775, 28759, 29323,
               <br/>
               30179, 29988, 30786, 31065, 30824, 30296, 29756, 29692, 29273,
               <br/>
               30276, 31432, 31365, 30612, 30549, 30394, 29769, 29893, 29167,
               <br/>
               27749, 27160, 26525, 28931, 28723, 27428, 27023, 26830, 28560,
               <br/>
               28481, 26174, 27534, 27807, 28102, 28162, 29195, 29445, 29414,
               <br/>
               30530, 30639, 30936, 31190, 31567, 31401, 29963, 30065, 30557,
               <br/>
               29746, 32022, 32144, 32364, 31245, 31310, 31194, 31073, 30297,
               <br/>
               29964, 28794, 28299, 29199, 28787, 28587, 26877, 28088, 26945,
               <br/>
               28171, 28221, 26107, 27486, 28269, 28237, 28909, 29137, 29579,
               <br/>
               30506, 30984, 31508, 31142, 31903, 32477, 31824, 30874, 30635,
               <br/>
               30995, 31228, 32296, 33258, 33124, 32020, 32558, 32372, 31312,
               <br/>
               30865, 30547, 28993, 29482, 28588, 28511, 28489, 28298, 27865,
               <br/>
               26416, 27777, 28404, 26173, 27675, 27924, 28942, 28988, 29395,
               <br/>
               30841, 30428, 31647, 31779, 31346, 32504, 31941, 31512, 31522,
               <br/>
               30967, 32097, 31463, 33465, 33772, 34025, 33199, 33584, 32971,
               <br/>
               32102, 31852, 31567, 28251, 28600, 27969, 27959, 28926, 27367,
               <br/>
               26470, 25399, 27132, 27473, 26302, 26840, 28622, 28693, 28884,
               <br/>
               30155, 30783, 30808, 31783, 31891, 31807, 31880, 32223, 31916,
               <br/>
               31367, 32534, 32578, 32222, 34255, 34489, 34933, 33944, 33727,
               <br/>
               33744, 32924, 32623, 27059, 26999, 27311, 27146, 27870, 28068,
               <br/>
               26486, 25557, 25274, 26854, 27404, 25477, 27672, 27779, 28609,
               <br/>
               29181, 29583, 30881, 30706, 31228, 31907, 31441, 32270, 32418,
               <br/>
               32063, 32689, 32015, 33182, 32615, 34646, 34807, 35701, 34530,
               <br/>
               34456, 34317, 33162, 27005, 26158, 26600, 27243, 27558, 27358,
               <br/>
               26672, 25497, 25364, 24445, 26322, 26557, 26246, 26980, 27574,
               <br/>
               28619, 28916, 30076, 30504, 31002, 31656, 31086, 30981, 31732,
               <br/>
               32197, 33410, 32661, 32934, 33732, 32989, 35102, 36020, 35746,
               <br/>
               34717, 35008, 34557, 26578, 25791, 25025, 25992, 27322, 26744,
               <br/>
               25933, 25600, 24882, 24607, 23851, 25414, 26726, 25848, 26657,
               <br/>
               28090, 28430, 29237, 29512, 30342, 30788, 30742, 30887, 30680,
               <br/>
               31663, 32947, 32720, 32569, 33003, 33661, 33385, 35604, 36353,
               <br/>
               36292, 35524, 35445, 25102, 25908, 24657, 24155, 25498, 26300,
               <br/>
               25319, 24712, 25185, 24206, 24252, 23300, 25523, 26077, 24676,
               <br/>
               26846, 27089, 28348, 28652, 29024, 30177, 29430, 30227, 30490,
               <br/>
               30251, 32676, 32829, 33249, 32917, 32749, 34048, 33930, 35604,
               <br/>
               35983, 36436, 35528, 21853, 22728, 23464, 23647, 23274, 24228,
               <br/>
               24142, 23150, 23234, 23570, 23463, 23008, 23044, 24673, 24815,
               <br/>
               24571, 26114, 26795, 27453, 27918, 28566, 29028, 29036, 29408,
               <br/>
               29708, 31281, 31694, 32249, 32863, 32920, 32813, 34429, 33516,
               <br/>
               35451, 36568, 36401, 21887, 21186, 21904, 22953, 22908, 22653,
               <br/>
               22691, 23057, 22495, 22564, 23185, 21999, 22761, 22002, 23560,
               <br/>
               24692, 23817, 25734, 25873, 26385, 27471, 27227, 27852, 28059,
               <br/>
               28999, 30199, 30203, 31559, 31969, 32431, 32623, 33404, 34053,
               <br/>
               33552, 35263, 35917, 20652, 20409, 19773, 20591, 21938, 21874,
               <br/>
               20768, 21388, 22076, 21677, 21642, 21690, 21887, 21756, 21085,
               <br/>
               23345, 23388, 23563, 24613, 24670, 25691, 25793, 26210, 27470,
               <br/>
               26891, 28839, 29182, 29795, 30914, 31614, 32709, 32623, 32350,
               <br/>
               33171, 32963, 34900, 18930, 19299, 18659, 18512, 20116, 20853,
               <br/>
               20265, 18989, 20147, 20370, 20298, 20107, 21342, 20514, 20197,
               <br/>
               20421, 21889, 22783, 21956, 23145, 24154, 24102, 24208, 25255,
               <br/>
               25791, 26887, 27584, 28184, 29219, 29867, 30593, 31731, 31369,
               <br/>
               31274, 32649, 32232, 16886, 16786, 17099, 17377, 16948, 18425,
               <br/>
               18753, 17764, 17691, 18147, 18951, 18429, 19602, 19385, 18819,
               <br/>
               19359, 18960, 20787, 21086, 20798, 22448, 22014, 22295, 23319,
               <br/>
               23785, 25594, 25771, 26441, 27255, 27913, 28464, 30524, 30401,
               <br/>
               30156, 30689, 31099, 14720, 14970, 15225, 15747, 15942, 15700,
               <br/>
               16399, 16338, 16683, 16084, 16968, 17272, 18082, 17987, 18108,
               <br/>
               18007, 18018, 17879, 19177, 19778, 19711, 20124, 20698, 21133,
               <br/>
               21433, 23137, 24013, 24647, 25198, 25871, 26412, 27775, 28512,
               <br/>
               28875, 29467, 29184, 12611, 12721, 12852, 13739, 13922, 14281,
               <br/>
               13624, 14162, 15104, 14786, 14803, 15171, 16447, 16020, 15893,
               <br/>
               17084, 16736, 16713, 16362, 17885, 18411, 17731, 18441, 19037,
               <br/>
               19394, 20886, 21221, 22383, 22521, 23524, 23921, 25418, 26109,
               <br/>
               26513, 27256, 27213, 10883, 10814, 10899, 11511, 11804, 12330,
               <br/>
               12460, 11658, 12770, 13134, 13352, 12847, 14159, 14327, 14205,
               <br/>
               14686, 15101, 15536, 15149, 14896, 16242, 16235, 15950, 16869,
               <br/>
               17504, 18602, 18927, 19747, 20593, 20950, 21536, 22655, 23161,
               <br/>
               23743, 24905, 25248, 8754, 9220, 9190, 9744, 10084, 10593, 10611,
               <br/>
               10355, 10692, 10918, 11933, 11496, 11819, 12497, 12359, 13029,
               <br/>
               13274, 13811, 13371, 13246, 13600, 14196, 14402, 14425, 15349,
               <br/>
               16633, 16630, 17336, 17715, 18741, 19364, 20159, 20694, 21118,
               <br/>
               22320, 22433), .Dim = as.integer(c(36, 30)))
               <br/>
               wa = structure(c(1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
               <br/>
               1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0, 0, 0, 0, 1, 1, 1, 1, 1, 1, 1,
               <br/>
               1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
               <br/>
               1, 1, 1, 1, 1, 0, 0, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
               <br/>
               1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
               <br/>
               0, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
               <br/>
               1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0, 1, 1, 1, 1,
               <br/>
               1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
               <br/>
               1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
               <br/>
               1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
               <br/>
               1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
               <br/>
               1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
               <br/>
               1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
               <br/>
               1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
               <br/>
               1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
               <br/>
               1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
               <br/>
               1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
               <br/>
               1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
               <br/>
               1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
               <br/>
               1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
               <br/>
               1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
               <br/>
               1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
               <br/>
               1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
               <br/>
               1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
               <br/>
               1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
               <br/>
               1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
               <br/>
               1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
               <br/>
               1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
               <br/>
               1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
               <br/>
               1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
               <br/>
               1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
               <br/>
               1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
               <br/>
               1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
               <br/>
               1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
               <br/>
               1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
               <br/>
               1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0,
               <br/>
               1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
               <br/>
               1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0, 0, 1, 1, 1, 1, 1,
               <br/>
               1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
               <br/>
               1, 1, 1, 1, 1, 1, 1, 1, 0, 0, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
               <br/>
               1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
               <br/>
               1, 1, 0, 0, 0, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
               <br/>
               1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0, 0, 0, 0,
               <br/>
               0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
               <br/>
               1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
               <br/>
               0, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
               <br/>
               1, 1, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 1,
               <br/>
               1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0,
               <br/>
               0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
               <br/>
               1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0,
               <br/>
               0, 0, 0, 0, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
               <br/>
               1, 1, 1, 1, 1, 1, 1, 1))
               <br/>
               aa&lt;-fit701(x,y,etx,dtx,wa)
               <br/>
               可是得到的结果，并不是fit701中list（）里面的参数。
              </p>
              <p>
               原始数据是这样排列的，可是我不懂如何调用？
              </p>
              <p>
               “qdata.usa” &lt;-
               <br/>
               structure(list(x = as.integer(c(60, 61, 62, 63, 64, 65, 66, 67,
               <br/>
               68, 69, 70, 71, 72, 73, 74, 75, 76, 77, 78, 79, 80, 81, 82, 83,
               <br/>
               84, 85, 86, 87, 88, 89)), y = as.integer(c(1968, 1969, 1970,
               <br/>
               1971, 1972, 1973, 1974, 1975, 1976, 1977, 1978, 1979, 1980, 1981,
               <br/>
               1982, 1983, 1984, 1985, 1986, 1987, 1988, 1989, 1990, 1991, 1992,
               <br/>
               1993, 1994, 1995, 1996, 1997, 1998, 1999, 2000, 2001, 2002, 2003
               <br/>
               )), etx = structure(c(856167, 866287, 903314, 875491, 880448,
               <br/>
               892288, 897373, 906129, 891513, 886044, 905800, 940496, 1021278,
               <br/>
               1062928, 1075375, 1037117, 1056018, 1052845, 1020308, 1024275,
               <br/>
               1035615, 1012013, 1011827, 994493, 949037, 939460, 933858, 995398,
               <br/>
               978947, 984941, 1037744, 1054359, 1119331, 1138132, 1215153,
               <br/>
               1367422, 813957, 832384, 841000, 886880, 865567, 876908, 896311,
               <br/>
               909091, 925398, 916932, 918861, 947455, 969148, 994121, 1035118,
               <br/>
               1047506, 1009979, 1028611, 1025401, 993776, 997940, 1009834,
               <br/>
               984232, 993307, 976781, 932267, 922196, 917577, 977659, 962298,
               <br/>
               968349, 1020164, 1066360, 1105452, 1123898, 1199917, 769980,
               <br/>
               793648, 809970, 825184, 871486, 850010, 861881, 882401, 895572,
               <br/>
               912479, 904679, 907038, 927024, 920673, 947893, 991340, 1007339,
               <br/>
               975370, 997039, 998294, 971480, 980136, 999158, 971062, 980626,
               <br/>
               964923, 920409, 911400, 906473, 966663, 952002, 958203, 1041131,
               <br/>
               1051307, 1090074, 1108149, 730887, 754131, 761967, 797520, 812424,
               <br/>
               859139, 837780, 850030, 870859, 884269, 901900, 893952, 900452,
               <br/>
               928482, 922297, 949057, 992613, 1008397, 976113, 997755, 998713,
               <br/>
               971933, 980666, 1003543, 975933, 986025, 969719, 925486, 915683,
               <br/>
               911581, 972711, 957676, 978987, 1025291, 1035410, 1073948, 694120,
               <br/>
               716735, 735620, 748868, 783125, 796414, 842754, 820393, 831159,
               <br/>
               851027, 863458, 880029, 879452, 900385, 928190, 922262, 948390,
               <br/>
               991946, 1007581, 975680, 997198, 998365, 971835, 981866, 1006136,
               <br/>
               979040, 989062, 973764, 928910, 919807, 916406, 977380, 957302,
               <br/>
               962735, 1008859, 1018693, 660800, 678637, 706601, 720541, 732391,
               <br/>
               765768, 778637, 824789, 801253, 811707, 831060, 843141, 862694,
               <br/>
               864906, 883432, 907341, 898438, 920579, 959249, 971271, 937508,
               <br/>
               955284, 959793, 939209, 949592, 974070, 947836, 958501, 943473,
               <br/>
               900478, 892028, 888601, 953199, 939883, 945434, 991349, 633685,
               <br/>
               643455, 666122, 688795, 701209, 711768, 744987, 757324, 803146,
               <br/>
               779215, 788633, 807738, 816708, 822067, 826637, 846701, 870571,
               <br/>
               864505, 886836, 926358, 939387, 908921, 928135, 932209, 912993,
               <br/>
               923548, 947407, 922890, 933141, 919059, 877065, 868563, 872243,
               <br/>
               934290, 921171, 927025, 610486, 614039, 629797, 645396, 666993,
               <br/>
               678635, 689230, 722362, 734266, 780124, 755703, 765708, 787301,
               <br/>
               796780, 803904, 809868, 831355, 856847, 852011, 875706, 916458,
               <br/>
               931259, 897058, 909180, 913663, 894941, 905034, 929745, 905727,
               <br/>
               916601, 902826, 861400, 868256, 853229, 914644, 901758, 583782,
               <br/>
               586915, 570692, 605096, 619979, 641386, 653640, 664590, 698043,
               <br/>
               710755, 757283, 733811, 743285, 752597, 759744, 763839, 766951,
               <br/>
               786160, 807024, 800260, 819592, 855612, 872167, 847508, 859081,
               <br/>
               863602, 845901, 856738, 880054, 858040, 868975, 855817, 856238,
               <br/>
               847663, 832930, 893913, 553129, 557079, 565585, 544330, 577926,
               <br/>
               592750, 614323, 627086, 638650, 672465, 685353, 731561, 705800,
               <br/>
               719835, 731507, 740038, 745870, 752289, 773836, 796830, 792965,
               <br/>
               814900, 855420, 862899, 838560, 850050, 854363, 837706, 848156,
               <br/>
               872633, 851424, 862324, 844657, 833686, 825860, 811629, 521770,
               <br/>
               524742, 546846, 531810, 505837, 533771, 543806, 559832, 567126,
               <br/>
               573232, 599359, 606474, 657413, 670764, 684652, 695130, 702755,
               <br/>
               708533, 714563, 736428, 757770, 754756, 777624, 829298, 836574,
               <br/>
               813103, 823991, 829034, 812563, 823632, 848179, 827664, 847845,
               <br/>
               819982, 810102, 802793, 493069, 493613, 489228, 521292, 509924,
               <br/>
               487698, 521330, 536585, 558895, 571906, 582857, 617382, 615889,
               <br/>
               629391, 643271, 656880, 667056, 674891, 680522, 687546, 709690,
               <br/>
               732227, 727327, 751182, 801706, 809020, 785910, 797618, 802053,
               <br/>
               786636, 798045, 821777, 798523, 821496, 794293, 785218, 464965,
               <br/>
               462899, 449521, 460001, 489586, 476858, 454388, 486414, 500451,
               <br/>
               520873, 531577, 541572, 579817, 590133, 601112, 611638, 621567,
               <br/>
               628331, 632384, 635025, 638853, 658325, 688165, 711376, 734762,
               <br/>
               784702, 791570, 770146, 781352, 786237, 771921, 783023, 786504,
               <br/>
               770509, 794034, 767594, 435941, 432637, 428927, 421251, 432833,
               <br/>
               464833, 454759, 433625, 468772, 485424, 508535, 522305, 524287,
               <br/>
               536570, 549664, 562947, 575846, 588445, 597690, 604940, 610374,
               <br/>
               618475, 629835, 638876, 661118, 682842, 728947, 736425, 716127,
               <br/>
               727165, 732611, 719175, 756705, 756305, 741418, 765096, 408692,
               <br/>
               402218, 407575, 398877, 390529, 402041, 434255, 424518, 404333,
               <br/>
               439415, 455364, 478352, 494175, 506181, 518832, 531733, 544936,
               <br/>
               558247, 570767, 580312, 587293, 594013, 600377, 608298, 617027,
               <br/>
               639073, 659894, 705499, 712462, 693485, 705067, 710352, 715681,
               <br/>
               725129, 725091, 711218, 380807, 375963, 380169, 378506, 369129,
               <br/>
               360766, 372938, 405468, 395893, 376042, 410407, 426790, 450828,
               <br/>
               463599, 476265, 488639, 501616, 515586, 529234, 542437, 551994,
               <br/>
               560445, 570224, 577377, 585864, 594420, 615387, 636575, 680315,
               <br/>
               687559, 670185, 681390, 695365, 682826, 692940, 693424, 348190,
               <br/>
               348880, 345547, 353216, 352107, 343718, 337013, 350924, 386492,
               <br/>
               378492, 356556, 393999, 407457, 419043, 431096, 442350, 452940,
               <br/>
               465081, 477747, 490662, 502117, 511231, 521203, 540528, 547729,
               <br/>
               556190, 564112, 585037, 604955, 647256, 655104, 638556, 651132,
               <br/>
               660748, 649399, 659990, 312353, 317343, 320291, 319364, 326669,
               <br/>
               325651, 318197, 312252, 326384, 362430, 354026, 333083, 368217,
               <br/>
               378617, 390303, 401778, 412498, 422828, 434639, 447633, 460160,
               <br/>
               472377, 481592, 491871, 510544, 517581, 525378, 534002, 553631,
               <br/>
               573151, 614325, 621876, 606071, 615676, 625752, 615365, 278317,
               <br/>
               282976, 281045, 297607, 297834, 307266, 308946, 303543, 298591,
               <br/>
               315332, 356824, 350445, 333174, 340792, 350879, 361371, 371677,
               <br/>
               381554, 390726, 402038, 413924, 426669, 441344, 462018, 472525,
               <br/>
               490868, 497465, 506050, 514196, 533545, 553352, 593222, 573575,
               <br/>
               570641, 580495, 590742, 246589, 250690, 245642, 257870, 272812,
               <br/>
               271793, 280268, 281212, 274785, 269800, 285253, 321181, 302394,
               <br/>
               308751, 315784, 324947, 334281, 343887, 353005, 361697, 371991,
               <br/>
               384214, 397314, 413871, 433410, 443319, 460297, 467521, 475505,
               <br/>
               483671, 503005, 521959, 528796, 537217, 535475, 545338, 218385,
               <br/>
               222159, 235718, 221534, 229560, 240281, 236160, 241397, 238536,
               <br/>
               229347, 221123, 232063, 262589, 268088, 274795, 280870, 289391,
               <br/>
               298159, 307301, 316436, 324212, 335173, 344376, 352360, 367467,
               <br/>
               384802, 393710, 409973, 416410, 423822, 431991, 449167, 476256,
               <br/>
               491351, 500206, 499428, 191831, 194665, 208358, 215780, 203574,
               <br/>
               214408, 229348, 228547, 237745, 238281, 230671, 224253, 225823,
               <br/>
               231064, 237512, 244519, 250579, 259603, 268728, 278865, 288280,
               <br/>
               297598, 307532, 315275, 322913, 336912, 352664, 361854, 376721,
               <br/>
               383001, 390718, 398318, 404709, 440534, 454061, 463344, 167849,
               <br/>
               170176, 170352, 188680, 194666, 180977, 191234, 205654, 203857,
               <br/>
               211889, 211198, 203570, 198428, 203255, 208777, 214984, 221582,
               <br/>
               227259, 236022, 245372, 254897, 265010, 273384, 279758, 286839,
               <br/>
               294185, 306825, 322103, 330490, 344466, 351051, 358174, 373147,
               <br/>
               369929, 405690, 417949, 145946, 148393, 143516, 151984, 169036,
               <br/>
               173413, 158822, 168511, 181773, 179350, 185783, 184670, 180549,
               <br/>
               183544, 188156, 192055, 196646, 201516, 205290, 212744, 219811,
               <br/>
               228050, 237226, 249171, 255501, 262325, 269022, 281504, 295494,
               <br/>
               303525, 317298, 323430, 319547, 340044, 336089, 371368, 125564,
               <br/>
               128851, 125396, 126958, 134037, 149410, 152316, 137075, 145796,
               <br/>
               157832, 154337, 159741, 158816, 161550, 165031, 169425, 172718,
               <br/>
               176646, 180713, 184009, 190366, 197346, 203495, 211185, 221884,
               <br/>
               227864, 233974, 240702, 251805, 264712, 272624, 285042, 280589,
               <br/>
               288454, 308106, 303998, 118294.174369962, 117593.628507295, 119786.337057443,
               <br/>
               121733.854555658, 118728.512804816, 129075.575196409, 131373.365625957,
               <br/>
               124444.967044179, 123933.568564432, 127128.057698194, 132760.446434037,
               <br/>
               129103.597030915, 137321, 134699, 138443, 142392, 146985, 150549,
               <br/>
               154402, 158689, 161412, 168523, 173121, 171462, 177988, 187182,
               <br/>
               192148, 197978, 203740, 213411, 225037, 231900, 247227, 250442,
               <br/>
               258129, 276966, 96295.436345537, 97930.8887291229, 99599.0501603805,
               <br/>
               103013.874737308, 104289.527596505, 102706.409689194, 107279.134553700,
               <br/>
               106880.084172105, 109137.008461453, 105218.464550382, 111001.424178741,
               <br/>
               112990.134277182, 118289, 117771, 115720, 118882, 122214, 126414,
               <br/>
               129315, 132401, 135164, 137533, 144002, 151735, 150320, 156241,
               <br/>
               164259, 169283, 174439, 179759, 188933, 199339, 209266, 218635,
               <br/>
               221409, 228887, 75689.0029792667, 76349.2036237612, 77135.4425731136,
               <br/>
               82459.0604973551, 83557.3942968323, 85712.0491275005, 81768.8507326564,
               <br/>
               84997.8320666383, 90651.550313127, 88742.970268134, 88845.0012768286,
               <br/>
               91053.6725238645, 98712, 100443, 100221, 98062, 100740, 103982,
               <br/>
               108035, 110707, 112906, 115577, 118031, 124972, 131717, 130542,
               <br/>
               135678, 143252, 147619, 152359, 157609, 165712, 174659, 183063,
               <br/>
               191560, 194338, 61829.1830637757, 61437.1759305036, 61920.0832685924,
               <br/>
               65397.0161028321, 67061.6261035384, 70049.9703368882, 70788.534501024,
               <br/>
               66232.1617345858, 72549.7259693481, 74617.7056289286, 75856.2209195565,
               <br/>
               72987.1832050286, 80441, 82144, 83849, 82960, 80274, 82298, 85108,
               <br/>
               88823, 90412, 92164, 93964, 101263, 107237, 113121, 112076, 117038,
               <br/>
               123569, 127522, 132180, 136825, 142566, 150887, 158556, 166688,
               <br/>
               47640.0101531432, 50176.0216600389, 50012.7591166765, 53027.6740841019,
               <br/>
               54877.9829088755, 57648.0040612573, 57745.9615872747, 56352.7878839157,
               <br/>
               58186.7704543532, 59416.6816143498, 64940.397664777, 62562.206616465,
               <br/>
               64320, 66906, 69284, 71351, 70635, 68463, 71073, 74925, 78939,
               <br/>
               81234, 82499, 79650, 85828, 90988, 95962, 95543, 99802, 105534,
               <br/>
               109371, 113474, 114786, 121844, 128943, 136238), .Dim = as.integer(c(36,
               <br/>
               30))), dtx = structure(c(21137, 20362, 21053, 20598, 20766, 21026,
               <br/>
               20378, 19689, 19332, 18437, 18790, 18119, 19212, 19062, 18826,
               <br/>
               18446, 18507, 18680, 17681, 17509, 17440, 16210, 15930, 15402,
               <br/>
               14878, 14538, 14379, 14708, 14141, 13954, 13893, 14016, 14152,
               <br/>
               14542, 15373, 16717, 21093, 21576, 20997, 21426, 22068, 21349,
               <br/>
               21138, 20328, 20488, 19787, 19608, 19382, 19207, 19856, 20094,
               <br/>
               19723, 19317, 19526, 19084, 18710, 18359, 17818, 16884, 16950,
               <br/>
               15977, 15912, 15208, 15199, 15354, 14830, 14780, 14881, 14657,
               <br/>
               14909, 15326, 16369, 22667, 22628, 22831, 22348, 23448, 23410,
               <br/>
               22314, 22219, 22310, 21906, 21298, 20757, 21006, 20407, 21555,
               <br/>
               21573, 21504, 21027, 21023, 20583, 19930, 19423, 19275, 18131,
               <br/>
               17868, 17670, 17082, 16338, 16482, 16170, 15974, 15506, 15832,
               <br/>
               15801, 16305, 16600, 23364, 22818, 23381, 23675, 23894, 24137,
               <br/>
               23362, 22804, 23007, 22916, 22873, 21796, 21637, 21943, 20956,
               <br/>
               22485, 22742, 22476, 21750, 21785, 21673, 20919, 20442, 19971,
               <br/>
               19180, 19154, 18385, 17748, 17441, 16617, 17078, 16626, 16511,
               <br/>
               16672, 16861, 16901, 23909, 23759, 23940, 23855, 25032, 24501,
               <br/>
               24242, 23920, 23576, 23459, 23996, 23420, 23125, 22804, 22615,
               <br/>
               22451, 23349, 23779, 23595, 22722, 22772, 22258, 21588, 21640,
               <br/>
               20894, 20503, 20443, 19449, 18778, 18121, 17514, 18254, 17744,
               <br/>
               17233, 17668, 17686, 25155, 24952, 25534, 25315, 25402, 26151,
               <br/>
               25495, 25265, 25308, 24760, 24643, 24737, 24891, 24157, 23878,
               <br/>
               24163, 23624, 24991, 25225, 24798, 24225, 23691, 23626, 22657,
               <br/>
               22407, 22503, 21216, 20935, 20201, 19320, 18908, 18538, 18614,
               <br/>
               18366, 18106, 18427, 23831, 24553, 24964, 25200, 25584, 25850,
               <br/>
               25891, 25057, 25424, 25273, 25329, 24966, 25689, 25347, 24881,
               <br/>
               24465, 24992, 24617, 25119, 25628, 25696, 24386, 24367, 24049,
               <br/>
               23170, 23352, 23163, 22152, 21759, 21043, 20462, 19743, 18960,
               <br/>
               19340, 18993, 18805, 27003, 24676, 25690, 25921, 26864, 26654,
               <br/>
               26479, 26549, 26216, 26176, 26621, 25567, 26743, 26416, 26337,
               <br/>
               26011, 25668, 26082, 25560, 26453, 26724, 26240, 25331, 25501,
               <br/>
               25109, 24941, 24594, 24254, 23058, 23098, 22022, 21115, 20658,
               <br/>
               19994, 20463, 20079, 27021, 26871, 25138, 25785, 26792, 26957,
               <br/>
               26819, 26882, 26967, 26236, 27039, 26615, 26988, 26943, 27258,
               <br/>
               27080, 26638, 26587, 26986, 26482, 27667, 27512, 26971, 25991,
               <br/>
               26531, 26854, 25485, 25253, 25091, 24331, 23777, 23073, 22082,
               <br/>
               21299, 20769, 21144, 25507, 27135, 27237, 25408, 26978, 27459,
               <br/>
               27157, 26608, 27436, 27428, 27163, 27181, 28148, 27404, 27783,
               <br/>
               28498, 28289, 27988, 27723, 27761, 27619, 28018, 28061, 28183,
               <br/>
               27011, 27795, 27550, 26821, 26551, 25998, 25024, 25060, 23621,
               <br/>
               23111, 22259, 21860, 27167, 26198, 27741, 27687, 26400, 27716,
               <br/>
               27258, 26867, 27577, 27489, 28053, 27411, 28736, 28704, 28172,
               <br/>
               29089, 28964, 29432, 28934, 28155, 28648, 27956, 28751, 29132,
               <br/>
               28810, 28816, 28703, 28552, 28017, 27833, 26962, 25849, 25690,
               <br/>
               24481, 24091, 23025, 27269, 26523, 25626, 27679, 28776, 26401,
               <br/>
               26874, 27039, 26966, 27466, 27921, 27959, 28297, 28977, 29096,
               <br/>
               29153, 29592, 30090, 30102, 29486, 29148, 29060, 28329, 29878,
               <br/>
               30110, 30663, 29723, 29800, 29578, 28528, 28667, 28204, 27007,
               <br/>
               26412, 25839, 24264, 29853, 27868, 27050, 26600, 29170, 29366,
               <br/>
               27147, 27534, 27931, 27909, 28493, 28336, 29775, 28759, 29323,
               <br/>
               30179, 29988, 30786, 31065, 30824, 30296, 29756, 29692, 29273,
               <br/>
               30276, 31432, 31365, 30612, 30549, 30394, 29769, 29893, 29167,
               <br/>
               27749, 27160, 26525, 28931, 28723, 27428, 27023, 26830, 28560,
               <br/>
               28481, 26174, 27534, 27807, 28102, 28162, 29195, 29445, 29414,
               <br/>
               30530, 30639, 30936, 31190, 31567, 31401, 29963, 30065, 30557,
               <br/>
               29746, 32022, 32144, 32364, 31245, 31310, 31194, 31073, 30297,
               <br/>
               29964, 28794, 28299, 29199, 28787, 28587, 26877, 28088, 26945,
               <br/>
               28171, 28221, 26107, 27486, 28269, 28237, 28909, 29137, 29579,
               <br/>
               30506, 30984, 31508, 31142, 31903, 32477, 31824, 30874, 30635,
               <br/>
               30995, 31228, 32296, 33258, 33124, 32020, 32558, 32372, 31312,
               <br/>
               30865, 30547, 28993, 29482, 28588, 28511, 28489, 28298, 27865,
               <br/>
               26416, 27777, 28404, 26173, 27675, 27924, 28942, 28988, 29395,
               <br/>
               30841, 30428, 31647, 31779, 31346, 32504, 31941, 31512, 31522,
               <br/>
               30967, 32097, 31463, 33465, 33772, 34025, 33199, 33584, 32971,
               <br/>
               32102, 31852, 31567, 28251, 28600, 27969, 27959, 28926, 27367,
               <br/>
               26470, 25399, 27132, 27473, 26302, 26840, 28622, 28693, 28884,
               <br/>
               30155, 30783, 30808, 31783, 31891, 31807, 31880, 32223, 31916,
               <br/>
               31367, 32534, 32578, 32222, 34255, 34489, 34933, 33944, 33727,
               <br/>
               33744, 32924, 32623, 27059, 26999, 27311, 27146, 27870, 28068,
               <br/>
               26486, 25557, 25274, 26854, 27404, 25477, 27672, 27779, 28609,
               <br/>
               29181, 29583, 30881, 30706, 31228, 31907, 31441, 32270, 32418,
               <br/>
               32063, 32689, 32015, 33182, 32615, 34646, 34807, 35701, 34530,
               <br/>
               34456, 34317, 33162, 27005, 26158, 26600, 27243, 27558, 27358,
               <br/>
               26672, 25497, 25364, 24445, 26322, 26557, 26246, 26980, 27574,
               <br/>
               28619, 28916, 30076, 30504, 31002, 31656, 31086, 30981, 31732,
               <br/>
               32197, 33410, 32661, 32934, 33732, 32989, 35102, 36020, 35746,
               <br/>
               34717, 35008, 34557, 26578, 25791, 25025, 25992, 27322, 26744,
               <br/>
               25933, 25600, 24882, 24607, 23851, 25414, 26726, 25848, 26657,
               <br/>
               28090, 28430, 29237, 29512, 30342, 30788, 30742, 30887, 30680,
               <br/>
               31663, 32947, 32720, 32569, 33003, 33661, 33385, 35604, 36353,
               <br/>
               36292, 35524, 35445, 25102, 25908, 24657, 24155, 25498, 26300,
               <br/>
               25319, 24712, 25185, 24206, 24252, 23300, 25523, 26077, 24676,
               <br/>
               26846, 27089, 28348, 28652, 29024, 30177, 29430, 30227, 30490,
               <br/>
               30251, 32676, 32829, 33249, 32917, 32749, 34048, 33930, 35604,
               <br/>
               35983, 36436, 35528, 21853, 22728, 23464, 23647, 23274, 24228,
               <br/>
               24142, 23150, 23234, 23570, 23463, 23008, 23044, 24673, 24815,
               <br/>
               24571, 26114, 26795, 27453, 27918, 28566, 29028, 29036, 29408,
               <br/>
               29708, 31281, 31694, 32249, 32863, 32920, 32813, 34429, 33516,
               <br/>
               35451, 36568, 36401, 21887, 21186, 21904, 22953, 22908, 22653,
               <br/>
               22691, 23057, 22495, 22564, 23185, 21999, 22761, 22002, 23560,
               <br/>
               24692, 23817, 25734, 25873, 26385, 27471, 27227, 27852, 28059,
               <br/>
               28999, 30199, 30203, 31559, 31969, 32431, 32623, 33404, 34053,
               <br/>
               33552, 35263, 35917, 20652, 20409, 19773, 20591, 21938, 21874,
               <br/>
               20768, 21388, 22076, 21677, 21642, 21690, 21887, 21756, 21085,
               <br/>
               23345, 23388, 23563, 24613, 24670, 25691, 25793, 26210, 27470,
               <br/>
               26891, 28839, 29182, 29795, 30914, 31614, 32709, 32623, 32350,
               <br/>
               33171, 32963, 34900, 18930, 19299, 18659, 18512, 20116, 20853,
               <br/>
               20265, 18989, 20147, 20370, 20298, 20107, 21342, 20514, 20197,
               <br/>
               20421, 21889, 22783, 21956, 23145, 24154, 24102, 24208, 25255,
               <br/>
               25791, 26887, 27584, 28184, 29219, 29867, 30593, 31731, 31369,
               <br/>
               31274, 32649, 32232, 16886, 16786, 17099, 17377, 16948, 18425,
               <br/>
               18753, 17764, 17691, 18147, 18951, 18429, 19602, 19385, 18819,
               <br/>
               19359, 18960, 20787, 21086, 20798, 22448, 22014, 22295, 23319,
               <br/>
               23785, 25594, 25771, 26441, 27255, 27913, 28464, 30524, 30401,
               <br/>
               30156, 30689, 31099, 14720, 14970, 15225, 15747, 15942, 15700,
               <br/>
               16399, 16338, 16683, 16084, 16968, 17272, 18082, 17987, 18108,
               <br/>
               18007, 18018, 17879, 19177, 19778, 19711, 20124, 20698, 21133,
               <br/>
               21433, 23137, 24013, 24647, 25198, 25871, 26412, 27775, 28512,
               <br/>
               28875, 29467, 29184, 12611, 12721, 12852, 13739, 13922, 14281,
               <br/>
               13624, 14162, 15104, 14786, 14803, 15171, 16447, 16020, 15893,
               <br/>
               17084, 16736, 16713, 16362, 17885, 18411, 17731, 18441, 19037,
               <br/>
               19394, 20886, 21221, 22383, 22521, 23524, 23921, 25418, 26109,
               <br/>
               26513, 27256, 27213, 10883, 10814, 10899, 11511, 11804, 12330,
               <br/>
               12460, 11658, 12770, 13134, 13352, 12847, 14159, 14327, 14205,
               <br/>
               14686, 15101, 15536, 15149, 14896, 16242, 16235, 15950, 16869,
               <br/>
               17504, 18602, 18927, 19747, 20593, 20950, 21536, 22655, 23161,
               <br/>
               23743, 24905, 25248, 8754, 9220, 9190, 9744, 10084, 10593, 10611,
               <br/>
               10355, 10692, 10918, 11933, 11496, 11819, 12497, 12359, 13029,
               <br/>
               13274, 13811, 13371, 13246, 13600, 14196, 14402, 14425, 15349,
               <br/>
               16633, 16630, 17336, 17715, 18741, 19364, 20159, 20694, 21118,
               <br/>
               22320, 22433), .Dim = as.integer(c(36, 30))), wa = structure(c(1,
               <br/>
               1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
               <br/>
               1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0, 0, 0, 0, 1, 1, 1, 1, 1, 1, 1,
               <br/>
               1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
               <br/>
               1, 1, 1, 1, 1, 0, 0, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
               <br/>
               1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
               <br/>
               0, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
               <br/>
               1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0, 1, 1, 1, 1,
               <br/>
               1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
               <br/>
               1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
               <br/>
               1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
               <br/>
               1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
               <br/>
               1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
               <br/>
               1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
               <br/>
               1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
               <br/>
               1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
               <br/>
               1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
               <br/>
               1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
               <br/>
               1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
               <br/>
               1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
               <br/>
               1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
               <br/>
               1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
               <br/>
               1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
               <br/>
               1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
               <br/>
               1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
               <br/>
               1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
               <br/>
               1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
               <br/>
               1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
               <br/>
               1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
               <br/>
               1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
               <br/>
               1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
               <br/>
               1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
               <br/>
               1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
               <br/>
               1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
               <br/>
               1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
               <br/>
               1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
               <br/>
               1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0,
               <br/>
               1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
               <br/>
               1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0, 0, 1, 1, 1, 1, 1,
               <br/>
               1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
               <br/>
               1, 1, 1, 1, 1, 1, 1, 1, 0, 0, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
               <br/>
               1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
               <br/>
               1, 1, 0, 0, 0, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
               <br/>
               1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0, 0, 0, 0,
               <br/>
               0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
               <br/>
               1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
               <br/>
               0, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
               <br/>
               1, 1, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 1,
               <br/>
               1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0,
               <br/>
               0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
               <br/>
               1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0,
               <br/>
               0, 0, 0, 0, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
               <br/>
               1, 1, 1, 1, 1, 1, 1, 1),
               <br/>
               .Dim = as.integer(c(36, 30)))),
               <br/>
               .Names = c(“x”, “y”, “etx”, “dtx”, “wa”))
              </p>
             </div>
             <!-- .bbp-reply-content -->
            </div>
            <!-- .reply -->
            <div class="bbp-reply-header" id="post-334550">
             <div class="bbp-meta">
              <span class="bbp-reply-post-date">
               2012年7月13日 下午12:00
              </span>
              <a class="bbp-reply-permalink" href="http://cos.name/cn/topic/107480/#post-334550">
               2 楼
              </a>
              <span class="bbp-admin-links">
              </span>
             </div>
             <!-- .bbp-meta -->
            </div>
            <!-- #post-334550 -->
            <div class="even bbp-parent-forum-999 bbp-parent-topic-107480 bbp-reply-position-2 user-id-109572 topic-author post-334550 reply type-reply status-publish hentry">
             <div class="bbp-reply-author">
              <a class="bbp-author-avatar" href="http://cos.name/cn/profile/109572/" rel="nofollow" title="查看ywh19860616的档案">
               <img src="http://sdn.geekzu.org/avatar/8599115e910c111374c57034b9f33bda?s=80&amp;d=monsterid&amp;r=g"/>
              </a>
              <br/>
              <a class="bbp-author-name" href="http://cos.name/cn/profile/109572/" rel="nofollow" title="查看ywh19860616的档案">
               ywh19860616
              </a>
              <br/>
              <div class="bbp-author-role">
               普通会员
              </div>
             </div>
             <!-- .bbp-reply-author -->
             <div class="bbp-reply-content">
              <p>
               谢谢各位了
               <br/>
               刚记起来好像最后用print命令可以输出结果。
              </p>
             </div>
             <!-- .bbp-reply-content -->
            </div>
            <!-- .reply -->
            <div class="bbp-reply-header" id="post-334553">
             <div class="bbp-meta">
              <span class="bbp-reply-post-date">
               2012年7月13日 下午12:28
              </span>
              <a class="bbp-reply-permalink" href="http://cos.name/cn/topic/107480/#post-334553">
               3 楼
              </a>
              <span class="bbp-admin-links">
              </span>
             </div>
             <!-- .bbp-meta -->
            </div>
            <!-- #post-334553 -->
            <div class="odd bbp-parent-forum-999 bbp-parent-topic-107480 bbp-reply-position-3 user-id-107675 post-334553 reply type-reply status-publish hentry">
             <div class="bbp-reply-author">
              <a class="bbp-author-avatar" href="http://cos.name/cn/profile/107675/" rel="nofollow" title="查看肖楠的档案">
               <img src="http://sdn.geekzu.org/avatar/66aa115466aeedf71937ff80c00b7982?s=80&amp;d=monsterid&amp;r=g"/>
              </a>
              <br/>
              <a class="bbp-author-name" href="http://cos.name/cn/profile/107675/" rel="nofollow" title="查看肖楠的档案">
               肖楠
              </a>
              <br/>
              <div class="bbp-author-role">
               版主
              </div>
             </div>
             <!-- .bbp-reply-author -->
             <div class="bbp-reply-content">
              <p>
               看文档啊少年 p11-p14 详细介绍了每个参数的数据类型和输出结果
              </p>
              <p>
               <a href="http://www.jpmorgan.com/directdoc/LM_User_Guide_V2.pdf" rel="nofollow">
                http://www.jpmorgan.com/directdoc/LM_User_Guide_V2.pdf
               </a>
              </p>
              <p>
              </p>
              <pre class="highlight ">source('http://perso.univ-rennes1.fr/arthur.charpentier/R/fitModels.r')
source('http://perso.univ-rennes1.fr/arthur.charpentier/R/simModels.r')
source('http://perso.univ-rennes1.fr/arthur.charpentier/R/sampleData.r')
x   = qdata.usa$x
y   = qdata.usa$y
etx = qdata.usa$etx
dtx = qdata.usa$dtx
wa  = qdata.usa$wa
res = fit707(x, y, etx, dtx, wa)
names(res)
#  [1] "kappa1"  "kappa2"  "kappa3"  "beta2"   "beta3"   "beta4v"  "gamma4"  "x0"
#  [9] "x"       "y"       "cy"      "epsilon" "mhat"    "mtx"     "npar"    "ll"
# [17] "BIC"     "wa"
</pre>
              <p>
              </p>
             </div>
             <!-- .bbp-reply-content -->
            </div>
            <!-- .reply -->
            <div class="bbp-reply-header" id="post-334554">
             <div class="bbp-meta">
              <span class="bbp-reply-post-date">
               2012年7月13日 下午12:37
              </span>
              <a class="bbp-reply-permalink" href="http://cos.name/cn/topic/107480/#post-334554">
               4 楼
              </a>
              <span class="bbp-admin-links">
              </span>
             </div>
             <!-- .bbp-meta -->
            </div>
            <!-- #post-334554 -->
            <div class="even bbp-parent-forum-999 bbp-parent-topic-107480 bbp-reply-position-4 user-id-109572 topic-author post-334554 reply type-reply status-publish hentry">
             <div class="bbp-reply-author">
              <a class="bbp-author-avatar" href="http://cos.name/cn/profile/109572/" rel="nofollow" title="查看ywh19860616的档案">
               <img src="http://sdn.geekzu.org/avatar/8599115e910c111374c57034b9f33bda?s=80&amp;d=monsterid&amp;r=g"/>
              </a>
              <br/>
              <a class="bbp-author-name" href="http://cos.name/cn/profile/109572/" rel="nofollow" title="查看ywh19860616的档案">
               ywh19860616
              </a>
              <br/>
              <div class="bbp-author-role">
               普通会员
              </div>
             </div>
             <!-- .bbp-reply-author -->
             <div class="bbp-reply-content">
              <p>
               回复 第3楼 的 nan.xiao：
               <br/>
               非常感谢版主，您这个太方便了。
               <br/>
               也非常感谢您提供的那个文章，刚才没有看到。
              </p>
             </div>
             <!-- .bbp-reply-content -->
            </div>
            <!-- .reply -->
            <div class="bbp-reply-header" id="post-334604">
             <div class="bbp-meta">
              <span class="bbp-reply-post-date">
               2012年7月15日 上午1:40
              </span>
              <a class="bbp-reply-permalink" href="http://cos.name/cn/topic/107480/#post-334604">
               5 楼
              </a>
              <span class="bbp-admin-links">
              </span>
             </div>
             <!-- .bbp-meta -->
            </div>
            <!-- #post-334604 -->
            <div class="odd bbp-parent-forum-999 bbp-parent-topic-107480 bbp-reply-position-5 user-id-109572 topic-author post-334604 reply type-reply status-publish hentry">
             <div class="bbp-reply-author">
              <a class="bbp-author-avatar" href="http://cos.name/cn/profile/109572/" rel="nofollow" title="查看ywh19860616的档案">
               <img src="http://sdn.geekzu.org/avatar/8599115e910c111374c57034b9f33bda?s=80&amp;d=monsterid&amp;r=g"/>
              </a>
              <br/>
              <a class="bbp-author-name" href="http://cos.name/cn/profile/109572/" rel="nofollow" title="查看ywh19860616的档案">
               ywh19860616
              </a>
              <br/>
              <div class="bbp-author-role">
               普通会员
              </div>
             </div>
             <!-- .bbp-reply-author -->
             <div class="bbp-reply-content">
              <p>
               回复 第3楼 的 nan.xiao：
               <br/>
               nan老师，您好
               <br/>
               我用自己数据代入， 总是有问题，请您帮忙看看。
               <br/>
               source(“fitModels.r”)
               <br/>
               source(“fit701.r”)
               <br/>
               xv&lt;-as.integer(21:90)
               <br/>
               xv&lt;-data.frame(xv)
               <br/>
               yv&lt;-as.integer(1995:2009)
               <br/>
               yv&lt;-data.frame(yv)
               <br/>
               etx&lt;-read.csv(‘etx.csv’,header = FALSE)
               <br/>
               dtx&lt;-read.csv(‘dtx.csv’,header = FALSE)
               <br/>
               wa&lt;-read.csv(‘wa.csv’,header = FALSE)
               <br/>
               aa&lt;-fit701(xv,yv,etx,dtx,wa)
              </p>
              <p>
               error ：only defined for equally-sized data frames
              </p>
             </div>
             <!-- .bbp-reply-content -->
            </div>
            <!-- .reply -->
            <div class="bbp-reply-header" id="post-334605">
             <div class="bbp-meta">
              <span class="bbp-reply-post-date">
               2012年7月15日 上午1:41
              </span>
              <a class="bbp-reply-permalink" href="http://cos.name/cn/topic/107480/#post-334605">
               6 楼
              </a>
              <span class="bbp-admin-links">
              </span>
             </div>
             <!-- .bbp-meta -->
            </div>
            <!-- #post-334605 -->
            <div class="even bbp-parent-forum-999 bbp-parent-topic-107480 bbp-reply-position-6 user-id-109572 topic-author post-334605 reply type-reply status-publish hentry">
             <div class="bbp-reply-author">
              <a class="bbp-author-avatar" href="http://cos.name/cn/profile/109572/" rel="nofollow" title="查看ywh19860616的档案">
               <img src="http://sdn.geekzu.org/avatar/8599115e910c111374c57034b9f33bda?s=80&amp;d=monsterid&amp;r=g"/>
              </a>
              <br/>
              <a class="bbp-author-name" href="http://cos.name/cn/profile/109572/" rel="nofollow" title="查看ywh19860616的档案">
               ywh19860616
              </a>
              <br/>
              <div class="bbp-author-role">
               普通会员
              </div>
             </div>
             <!-- .bbp-reply-author -->
             <div class="bbp-reply-content">
              <p>
               回复 第3楼 的 nan.xiao：
              </p>
              <p>
               还有，在R2.14.2运行
               <br/>
               &gt; source(“fitModels.r”)
               <br/>
               警告信息：
               <br/>
               1: In grepl(“\n”, lines, fixed = TRUE) : 输入的字符串4不适用于此语言环境
               <br/>
               2: In grepl(“\n”, lines, fixed = TRUE) : 输入的字符串5不适用于此语言环境
               <br/>
               3: In grepl(“\n”, lines, fixed = TRUE) : 输入的字符串6不适用于此语言环境
               <br/>
               4: In grepl(“\n”, lines, fixed = TRUE) : 输入的字符串7不适用于此语言环境
               <br/>
               5: In grepl(“\n”, lines, fixed = TRUE) : 输入的字符串19不适用于此语言环境
               <br/>
               6: In grepl(“\n”, lines, fixed = TRUE) : 输入的字符串20不适用于此语言环境
               <br/>
               7: In grepl(“\n”, lines, fixed = TRUE) : 输入的字符串21不适用于此语言环境
               <br/>
               8: In grepl(“\n”, lines, fixed = TRUE) : 输入的字符串23不适用于此语言环境
               <br/>
               9: In grepl(“\n”, lines, fixed = TRUE) : 输入的字符串36不适用于此语言环境
               <br/>
               10: In grepl(“\n”, lines, fixed = TRUE) : 输入的字符串41不适用于此语言环境
               <br/>
               &gt;
               <br/>
               但是在R2.12.2就可以，为什么？
              </p>
             </div>
             <!-- .bbp-reply-content -->
            </div>
            <!-- .reply -->
            <div class="bbp-reply-header" id="post-334615">
             <div class="bbp-meta">
              <span class="bbp-reply-post-date">
               2012年7月15日 上午8:01
              </span>
              <a class="bbp-reply-permalink" href="http://cos.name/cn/topic/107480/#post-334615">
               7 楼
              </a>
              <span class="bbp-admin-links">
              </span>
             </div>
             <!-- .bbp-meta -->
            </div>
            <!-- #post-334615 -->
            <div class="odd bbp-parent-forum-999 bbp-parent-topic-107480 bbp-reply-position-7 user-id-107675 post-334615 reply type-reply status-publish hentry">
             <div class="bbp-reply-author">
              <a class="bbp-author-avatar" href="http://cos.name/cn/profile/107675/" rel="nofollow" title="查看肖楠的档案">
               <img src="http://sdn.geekzu.org/avatar/66aa115466aeedf71937ff80c00b7982?s=80&amp;d=monsterid&amp;r=g"/>
              </a>
              <br/>
              <a class="bbp-author-name" href="http://cos.name/cn/profile/107675/" rel="nofollow" title="查看肖楠的档案">
               肖楠
              </a>
              <br/>
              <div class="bbp-author-role">
               版主
              </div>
             </div>
             <!-- .bbp-reply-author -->
             <div class="bbp-reply-content">
              <p>
               <code>
                wa.csv
               </code>
               有问题，第92行，以及前面空白的若干行。
              </p>
             </div>
             <!-- .bbp-reply-content -->
            </div>
            <!-- .reply -->
            <div class="bbp-reply-header" id="post-334616">
             <div class="bbp-meta">
              <span class="bbp-reply-post-date">
               2012年7月15日 上午8:06
              </span>
              <a class="bbp-reply-permalink" href="http://cos.name/cn/topic/107480/#post-334616">
               8 楼
              </a>
              <span class="bbp-admin-links">
              </span>
             </div>
             <!-- .bbp-meta -->
            </div>
            <!-- #post-334616 -->
            <div class="even bbp-parent-forum-999 bbp-parent-topic-107480 bbp-reply-position-8 user-id-109572 topic-author post-334616 reply type-reply status-publish hentry">
             <div class="bbp-reply-author">
              <a class="bbp-author-avatar" href="http://cos.name/cn/profile/109572/" rel="nofollow" title="查看ywh19860616的档案">
               <img src="http://sdn.geekzu.org/avatar/8599115e910c111374c57034b9f33bda?s=80&amp;d=monsterid&amp;r=g"/>
              </a>
              <br/>
              <a class="bbp-author-name" href="http://cos.name/cn/profile/109572/" rel="nofollow" title="查看ywh19860616的档案">
               ywh19860616
              </a>
              <br/>
              <div class="bbp-author-role">
               普通会员
              </div>
             </div>
             <!-- .bbp-reply-author -->
             <div class="bbp-reply-content">
              <p>
               回复 第7楼 的 nan.xiao：
               <br/>
               非常抱歉，自己粗心了
               <br/>
               老师，第6楼的错误可以不用更换R版本解决吗？
              </p>
             </div>
             <!-- .bbp-reply-content -->
            </div>
            <!-- .reply -->
            <div class="bbp-reply-header" id="post-334617">
             <div class="bbp-meta">
              <span class="bbp-reply-post-date">
               2012年7月15日 上午8:29
              </span>
              <a class="bbp-reply-permalink" href="http://cos.name/cn/topic/107480/#post-334617">
               9 楼
              </a>
              <span class="bbp-admin-links">
              </span>
             </div>
             <!-- .bbp-meta -->
            </div>
            <!-- #post-334617 -->
            <div class="odd bbp-parent-forum-999 bbp-parent-topic-107480 bbp-reply-position-9 user-id-107675 post-334617 reply type-reply status-publish hentry">
             <div class="bbp-reply-author">
              <a class="bbp-author-avatar" href="http://cos.name/cn/profile/107675/" rel="nofollow" title="查看肖楠的档案">
               <img src="http://sdn.geekzu.org/avatar/66aa115466aeedf71937ff80c00b7982?s=80&amp;d=monsterid&amp;r=g"/>
              </a>
              <br/>
              <a class="bbp-author-name" href="http://cos.name/cn/profile/107675/" rel="nofollow" title="查看肖楠的档案">
               肖楠
              </a>
              <br/>
              <div class="bbp-author-role">
               版主
              </div>
             </div>
             <!-- .bbp-reply-author -->
             <div class="bbp-reply-content">
              <p>
               警告是他注释中的引号不是正常的英文引号导致的，把前方或后方那一大段注释的协议删了就行了。
              </p>
             </div>
             <!-- .bbp-reply-content -->
            </div>
            <!-- .reply -->
            <div class="bbp-reply-header" id="post-334618">
             <div class="bbp-meta">
              <span class="bbp-reply-post-date">
               2012年7月15日 上午8:40
              </span>
              <a class="bbp-reply-permalink" href="http://cos.name/cn/topic/107480/#post-334618">
               10 楼
              </a>
              <span class="bbp-admin-links">
              </span>
             </div>
             <!-- .bbp-meta -->
            </div>
            <!-- #post-334618 -->
            <div class="even bbp-parent-forum-999 bbp-parent-topic-107480 bbp-reply-position-10 user-id-109572 topic-author post-334618 reply type-reply status-publish hentry">
             <div class="bbp-reply-author">
              <a class="bbp-author-avatar" href="http://cos.name/cn/profile/109572/" rel="nofollow" title="查看ywh19860616的档案">
               <img src="http://sdn.geekzu.org/avatar/8599115e910c111374c57034b9f33bda?s=80&amp;d=monsterid&amp;r=g"/>
              </a>
              <br/>
              <a class="bbp-author-name" href="http://cos.name/cn/profile/109572/" rel="nofollow" title="查看ywh19860616的档案">
               ywh19860616
              </a>
              <br/>
              <div class="bbp-author-role">
               普通会员
              </div>
             </div>
             <!-- .bbp-reply-author -->
             <div class="bbp-reply-content">
              <p>
               回复 第9楼 的 nan.xiao：
              </p>
              <p>
               非常感谢老师，删除后没有那个错误提示了。
               <br/>
               谢谢。
              </p>
             </div>
             <!-- .bbp-reply-content -->
            </div>
            <!-- .reply -->
           </li>
           <!-- .bbp-body -->
           <li class="bbp-footer">
            <div class="bbp-reply-author">
             作者
            </div>
            <div class="bbp-reply-content">
             帖子
            </div>
            <!-- .bbp-reply-content -->
           </li>
           <!-- .bbp-footer -->
          </ul>
          <!-- #topic-107480-replies -->
          <div class="bbp-pagination">
           <div class="bbp-pagination-count">
            查看 10 个帖子 - 1 到 10（总计 10 个）
           </div>
           <div class="bbp-pagination-links">
           </div>
          </div>
          <div class="bbp-no-reply" id="no-reply-107480">
           <div class="bbp-template-notice">
            <p>
             您必须先登录才能回复该主题。
            </p>
           </div>
          </div>
         </div>
        </div>
        <!-- .entry-content -->
        <footer class="entry-meta">
        </footer>
        <!-- .entry-meta -->
       </article>
       <!-- #post -->
       <div class="comments-area" id="comments">
       </div>
       <!-- #comments -->
      </div>
      <!-- #content -->
     </div>
     <!-- #primary -->
     <div class="sidebar-container" id="tertiary" role="complementary">
      <div class="sidebar-inner">
       <div class="widget-area">
        <aside class="widget bbp_widget_login" id="bbp_login_widget-2">
         <h3 class="widget-title">
          登录
         </h3>
         <form action="http://cos.name/cn/wp-login.php" class="bbp-login-form" method="post">
          <fieldset>
           <legend>
            登录
           </legend>
           <div class="bbp-username">
            <label for="user_login">
             用户名:
            </label>
           </div>
           <div class="bbp-password">
            <label for="user_pass">
             密码:
            </label>
           </div>
           <div class="bbp-remember-me">
            <label for="rememberme">
             记住用户名
            </label>
           </div>
           <div class="bbp-submit-wrapper">
            <button class="button submit user-submit" id="user-submit" name="user-submit" tabindex="104" type="submit">
             登录
            </button>
           </div>
           <div class="bbp-login-links">
            <a class="bbp-register-link" href="http://cos.name/cn/wp-login.php?action=register" title="注册">
             注册
            </a>
            <a class="bbp-lostpass-link" href="http://cos.name/cn/wp-login.php?action=lostpassword" title="忘记密码">
             忘记密码
            </a>
           </div>
          </fieldset>
         </form>
        </aside>
        <aside class="widget widget_text" id="text-7">
         <h3 class="widget-title">
          搜索
         </h3>
         <div class="textwidget">
          <form action="http://www.google.com/search" id="bbp-search-form" method="get" onsubmit="Gsitesearch(this)" role="search">
           <div>
           </div>
          </form>
          <form id="bbp-search-form-baidu" onsubmit="g(this)" role="search">
           <div>
           </div>
          </form>
         </div>
        </aside>
        <aside class="widget widget_text" id="text-2">
         <h3 class="widget-title">
          新鲜事
         </h3>
         <div class="textwidget">
          <ul>
           <li>
            <a href="http://cos.name/cn/topics/">
             最新帖子
            </a>
           </li>
           <li>
            <a href="http://cos.name/cn/view/popular/">
             最热门主题
            </a>
           </li>
           <li>
            <a href="http://cos.name/cn/view/no-replies/">
             消灭零回复
            </a>
           </li>
          </ul>
         </div>
        </aside>
        <aside class="widget widget_text" id="text-3">
         <h3 class="widget-title">
          RSS订阅
         </h3>
         <div class="textwidget">
          <ul>
           <li>
            <img src="http://cos.name/wp-includes/images/rss.png"/>
            <a href="http://cos.name/cn/topics/feed/">
             所有主题
            </a>
           </li>
           <li>
            <img src="http://cos.name/wp-includes/images/rss.png"/>
            <a href="http://cos.name/cn/forums/feed/">
             所有帖子
            </a>
           </li>
          </ul>
         </div>
        </aside>
       </div>
       <!-- .widget-area -->
      </div>
      <!-- .sidebar-inner -->
     </div>
     <!-- #tertiary -->
    </div>
    <!-- #main -->
    <footer class="site-footer" id="colophon" role="contentinfo">
     <div class="site-info">
      版权所有 © 2014 统计之都 | 由
      <a href="http://wordpress.org/">
       WordPress
      </a>
      构建 | 主题修改自
      <a href="http://wordpress.org/themes/twentythirteen">
       Twenty Thirteen
      </a>
     </div>
     <!-- .site-info -->
    </footer>
    <!-- #colophon -->
   </div>
   <!-- #page -->
  </div>
 </body>
</html>